Quasi-Patterns in a Model of Multi-Resonantly Forced Chemical Oscillations 
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Multi-frequency forcing of systems undergoing a Hopf bifurcation to spatially homogeneous os- 
cillations is investigated using a complex Ginzburg-Landau equation that systematically captures 
weak forcing functions that simultaneously hit the 1:1-, the 1:2-, and the l:3-resonance. Weakly 
nonlinear analysis shows that generically the forcing function can be tuned such that resonant triad 
interactions with weakly damped modes stabilize subharmonic quasipatterns with 4-fold and 5- 
fold rotational symmetry. In simulations starting from random initial conditions domains of these 
quasi-patterns compete and yield complex, slowly ordering patterns. 
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In recent years complex, but ordered spatio- 
temporal patterns characterized by multiple length 
scales have found considerable interest. In partic- 
ular on the surface of vertically vibrated fluid lay- 
ers (Faraday system) various kinds of fascinating 
periodic superlattice patterns and quasi-patterns 
have been found experimentally [l[ . Subsequently, 
such patterns have also been observed in optical 
systems in vertically vibrated fluid convection 
[3j, and on the surface of ferrofluids driven by 
time-periodic magnetic fields [4|. Here we show 
that superlattices and quasipatterns should be ac- 
cessible quite generally in a different class of sys- 
tems: resonantly forced systems undergoing a Hopf 
bifurcation to spatially homogeneous oscillations. 
Paradigmatic for such systems are chemical oscil- 
lations [5] . In chemical systems patterns with mul- 
tiple length scales have so far been obtained only 
by imposing an external length through spatially 
periodic illumination 

The stability of multi-mode patterns depends on 
the interaction between their constitutive Fourier 
modes. For small angles 9 between the modes 
the cross-coupling coefficient b(9) is twice as large 
as the self-coupling coefficient bo. Nearly paral- 
lel modes therefore suppress each other and un- 
less the cross-coupling coefficient decreases sub- 
stantially with increasing 9 only stripe-like pat- 
terns are stable. Strong angle dependence can arise 
if the basic modes couple to weakly damped, res- 
onating modes 0, § ■ Complex patterns with dif- 
ferent symmetries can then become stable @. 

The resonances stabilizing complex patterns 
have been studied in great detail in the Faraday 
system. Their spatio-temporal nature allows a 
very controled tuning through the frequency con- 
tent of the driving Broadly speaking, there 
are two mechanisms by which complex patterns 
can be stabilized: either by enhancing the self- 
damping bo d, [H, [lU or by reducing the cross- 
coupling coefficient b(9) [Til. Il3j. 



In this paper we exploit spatio-temporal reso- 
nances to induce complex spatial patterns in two- 
dimensional systems undergoing a Hopf bifurcation 
to spatially homogeneous oscillations. To this end 
we apply spatially homogeneous, resonant multi- 
frequency forcing. By including in the forcing spec- 
trum a frequency component close to twice the 
Hopf frequency (l:2-resonance) we excite standing 
waves with a wavenumber determined by the de- 
tuning between the forcing and the Hopf frequency 
[1J| . A second frequency near three times the Hopf 
frequency (l:3-resonance) induces a quadratic in- 
teraction term, which otherwise is not allowed 
in the normal form for Hopf bifurcations. To 
avoid transcritical bifurcations off the basic state 
to hexagonal patterns we further include a second 
forcing frequency close to the l:2-resonance with 
a slightly different detuning. Within the weakly 
nonlinear regime we show that quite independent 
of the two parameters characterizing the unforced 
Hopf bifurcation the forcing function can be tuned 
such that instead of the usual stripe, spiral, and 
labyrinthine patterns [15[ one obtains superlattices 
and quasipatterns. 

The systematic, weakly nonlinear description of 
a weakly forced super-critical Hopf bifurcation is 
given by the complex Ginzburg-Landau equation 
for the complex oscillation amplitude A, which is 
extended to include near-resonant components of 
the forcing function 

fj A 

— = {l + if3)V 2 A+(Li + ia-{l + ia)\A\ 2 )A 

+-f(cosx + sm X e wt )A*+pe l ' s 'A* 2 . (1) 

Here \ measures the relative contributions from 
the two forcing components that are close to the 
l:2-resonance, which differ in their frequencies by 
v. The detuning between the Hopf frequency and 
half the frequency corresponding to the forcing 
7 cos x is given by a. The strength and phase of the 
l:3-forcing is given by p and respectively. Non- 
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linear interactions of the l:2-resonant forcing and 
the l:3-resonant forcing introduce an additional 
forcing near the Hopf frequency itself. To can- 
cel the resulting additional, inhomogeneous term 
in ([1]) we assume a further explicit forcing compo- 
nent near the l:l-resonance. It is straightforward 
to derive ([1]) from Oregonator-type models for the 
photosensitive Belousov-Zhabotinsky reaction 

The slight detuning between the two l:2-forcing 
components introduces the explicit, periodic time 
dependence of the coefficients in ([T]). Using Flo- 
quet theory we determine the instability of the ba- 
sic state A = with respect to time-periodic solu- 
tions that are phase- locked to the forcing [lTj . Due 
to the dispersion (3 the detuning a induces phase- 
locked modes with a non-zero wave number [LH ]. 
Depending on the forcing parameter \ and the de- 
tuning v the mode that destabilizes the basic state 
first is either harmonic or subharmonic relative to 
the period 2tt/v. 

A typical set of neutral curves ^ H ' SH ^{)i) for 
the harmonic and subharmonic mode is shown in 
FigfT] We focus here on the subharmonic case for 
which the time-shift symmetry t — > t + 2n/v sup- 
presses any transcritical bifurcations off A — to 
hexagons. Weakly damped harmonic modes ex- 
cited at quadratic order modify the competition 
between subharmonic modes of different orienta- 
tion [l(| E3] • Their effect is strongest if the forcing 
7 is only slightly below the critical forcing strength 
7c of the harmonic modes (inset of Fig[T]). For 
this reason we tune the forcing parameter \ an( A 
the detuning v so that 7c is only slightly above 
the critical forcing strength of the subhar- 

monic modes. In this paper we focus on the en- 
hancement of the self-damping bo of the subhar- 
monic modes and choose the forcing function such 

( H} 

that the critical wavenumber k c of the harmonic 
mode is close to twice that of the subharmonic 
mode, K = kc /kf 11 " 1 ~ 2. To reduce the compe- 
tition between modes subtending a specific angle 9 
a wave- number ratio K < 2 would be chosen [l7| • 
To compute the interaction between modes of 
different orientation within weakly nonlinear anal- 
ysis we expand the oscillation amplitude as (e ^ 1) 

A = e (A 1 e l ^ r + A 2 e l ^ r ) F(t) + h.o.t. (2) 

For subharmonic patterns F{t) has periodicity 
47r jv and the amplitude equations for A\p do not 
contain any quadratic terms, 



^- = XA 1 -b A 1 \A 1 \ 2 -b(e)A 1 \A 2 \ 2 , (3) 
with a similar equation for A 2 . 
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Figure 1: Neutral Curves for subharmonic (thick lines) 
and harmonic modes (thin) for p, = — 1, a = 4, yS = 3, 
X = 0.476718, and v = 4.2. 
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Figure 2: Angle-dependence of the effective cross- 
coupling b(8)/bo for a — —1 and $ = ^p. Other pa- 
rameters as in Fig[T] 



Relevant for the pattern selection is the ra- 
tio b(9)/bo. It is strongly affected by spatio- 
temporally resonant triads, which are induced by 
the l:3-forcing pe l *. The resulting p— dependence 
of b(9)/bo is shown in FigJ5J The stability con- 
ditions for rectangular patterns (corresponding to 
a rhombic arrangement of the wave vectors) are 
bo > and \b(9)/bo\ < 1. Thus, with increas- 
ing 1:3- forcing p a large range of angles arises 
for which rectangular patterns are stable, whereas 
without that forcing only stripe patterns would be 
obtained. 

Given b(9), the linear stability of various types of 
periodic super-lattice patterns comprised of three 
or more modes on a fixed periodic Fourier lattice 
can be determined systematically [llj]. The com- 
petition between different planforms involves, how- 
ever, often a collection of modes that cannot be 
represented on a single Fourier lattice. We have 
determined the relative stability of such planforms 
approximately by keeping all modes on the critical 
circle that are involved in the pattern competition 
and find bistability between a number of different 
complex patterns [17| . This approach ignores pos- 
sible side-band instabilities [19( and higher-order 
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Figure 3: Energy of planforms comprised of N equally- 
spaced Fourier modes. Parameters as in Figf2] Filled 
(open) symbols: linearly stable (unstable) planforms. 



resonances [20j and does not account for possible 
convergence problems due to small divisors [2l[ . 

To address the competition between different, si- 
multaneously stable planforms we exploit the vari- 
ational character of = — §§r for j = 

1, . . . , N. Fig. [3] shows the energies Fn of patterns 
comprised of N modes that are equally spaced on 
the critical circle as a function of the 1:3- forcing 
strength p. When starting from random initial 
conditions planforms with lower energy are ex- 
pected to invade those with higher energy. Thus, 
for p < 0.82 the final state is expected to consist 
of stripes, whereas for p > 1.07 patterns with four 
or more modes should dominate. 

To test the predictions of our weakly nonlinear 
analysis we have performed direct simulations of 
the complex Ginzburg-Landau equation ([T]). For 
small sizes they confirm the linear stability of pe- 
riodic patterns comprised of four modes. To in- 
vestigate the dependence of the pattern selection 
on the l:3-forcing strength p we performed sim- 
ulations in a large square system of linear size 
L m 473.39, which is equivalent to forty wave- 
lengths, for increasing forcing strengths p, start- 
ing from identical random initial conditions. We 
have chosen the system size such that the modes 
making up hexagons and super-squares have equal 
growth rates to bring out clearly how an increase 
in p alone tips the balance from hexagons to 4-fold 
patterns. For p = 0.9 a pattern with hexagonal 
planform rather than a stripe structure arises. Be- 
cause of the reflection symmetry A — > — A induced 
by the time-shift symmetry t — > t + lix jv domains 
with up- and down-hexagons coexist separated by 
walls containing narrow layers of triangle patterns 

Increasing p decreases the #-range over which 
modes suppress each other (cf. FiglJ) and more 
modes persist, as shown in FigQ] for p = 1.2. 
The pattern exhibits domains of periodic super- 




Figure 4: Partial view (size 0.5L x 0.35L) of four-mode 
pattern for p — 1.2. Other parameters as in Fig. [2] 




Figure 5: Partial view (size 0.5L x 0.35L) of five-mode 
pattern for p — 3. Other parameters as in Fig(2] 



lattice patterns with four-fold rotational symmetry 
('super-squares' and 'anti-squares' [HI], marked by 
dashed-dotted and dashed circles) as well as ele- 
ments with eight-fold rotational symmetry (solid 
circle). Increasing the l:3-forcing to p = 3 in- 
creases the number of persisting modes further and 
introduces numerous elements with five- and with 
ten-fold rotational symmetry (dashed and solid cir- 
cles) in Fig. [5] 

The patterns shown in Figs l4l5l are still evolv- 
ing, albeit very slowly. Nevertheless, it is clear 
that for p > 1.1 they will not evolve to simple 
hexagon states. While in our simulations for all 
values of p domains of hexagons appeared for early 
times, they were replaced for p > 1.1 by domains of 
patterns comprised of four or more modes, which 
have lower energy. A condensed view of the tem- 
poral evolution of the patterns for different val- 
ues of the forcing p is given in Fig [5] It depicts 
the evolution of the relevant number of Fourier 
modes e s , estimated by the spectral pattern en- 
tropy S = Pij^ n Pij- Here pij denotes the 
normalized power spectrum. Clearly the number 
of significant modes increases with p, albeit not 
monotonically at all times. 

In conclusion, we have shown that in systems 
undergoing a Hopf bifurcation to spatially homoge- 
neous oscillations multi-frequency forcing can sub- 
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Figure 6: Temporal evolution of the relevant number of 
Fourier modes e s for various values of the l:3-forcing 
strength p. Other parameters as in Fig[2] 

stantially reduce the competition between modes 
of different orientation leading to complex multi- 
mode patterns. By an appropriate choice of the 
amplitudes and phases of the forcing function, 
which constitute external control parameters, this 
regime should be accessible generically, essentially 
independent of the specifics of the unforced sys- 
tem. Our results should therefore apply to realis- 
tic chemical oscillators 0, [TH, [22| • From a practical 
point of view it should be mentioned, however, that 
the complex patterns possibly arise only very close 



to onset. This may require systems with relatively 
large aspect ratios and a very precise tuning of the 
forcing parameters. 

Using direct simulations of the complex 
Ginzburg-Landau equation we confirmed that 
these complex patterns arise from general ran- 
dom initial conditions. The appropriate, quanti- 
tative characterization of the transients, in which 
multi-mode structures like super-squares and anti- 
squares compete with each other, is still an open 
problem (cf. [Hj]). Another interesting question is 
the long-time scaling of the ordering of such com- 
plex structures. 

Compared to the Faraday system, the forced 
Hopf bifurcation considered here allows an addi- 
tional level of complexity by going slightly above 
the Hopf bifurcation. There complex patterns 
would compete with spatially homogeneous oscil- 
lations. For single-frequency forcing labyrinthine 
stripe patterns arise from the oscillations through 
front instabilities and stripe nucleation [22j |. It is 
not known what happens to this scenario when the 
stripes are unstable to the more complex patterns 
discussed here. 
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